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SOME APPLICATIONS OF WEIGHTED NORM INEQUALITIES TO 
THE ERROR ANALYSIS OF PDE CONSTRAINED OPTIMIZATION 

PROBLEMS* 

HARBIR ANTILt, ENRIQUE OTAROLAt, AND ABNER J. SALGADO§ 

Abstract. The purpose of this work is to illustrate how the theory of Muckenhoupt weights, 
Muckenhoupt weighted Sobolev spaces and the corresponding weighted norm inequalities can be used 
in the analysis and discretization of PDE constrained optimization problems. We consider: a linear 
quadratic constrained optimization problem where the state solves a nonuniformly elliptic equation; 
a problem where the cost involves pointwise observations of the state and one where the state has 
singular sources, e.g. point masses. For all three examples we propose and analyze numerical schemes 
and provide error estimates in two and three dimensions. While some of these problems might have 
been considered before in the literature, our approach allows for a simpler, Hilbert space-based, 
analysis and discretization and further generalizations. 
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1. Introduction. The purpose of this work is to show how the theory of Muck¬ 
enhoupt weights, Muckenhoupt weighted Sobolev spaces and weighted norm inequal¬ 
ities can be applied to analyze PDE constrained optimization problems, and their 
discretizations. These tools have already been shown to be essential in the analysis 
and discretization of problems constrained by equations involving fractional deriva¬ 
tives both in space and time 016], and here we extend their use to a new class of 
problems. 

We consider three illustrative examples. While some of them have been considered 
before, the techniques that we present are new and we believe they provide simpler 
arguments and allow for further generalizations. To describe them let D be an open 
and bounded polytopal domain of M" (n G {2, 3}) with Lipschitz boundary dVl. We 
will be concerned with the following problems: 

• Optimization with nonuniformly elliptic equations. Let w be a weight, that 
is, a positive and locally integrable function and yd G L^(w, D). Given a regulariza¬ 
tion parameter A > 0, we define the cost functional 

JAiy,u) yd\\‘L2(uj,n) + (^-l) 

We are then interested in finding min J 4 subject to the nonuniformly elliptic prob¬ 
lem 


— div(AVi/) = u in D, ?/ = 0 on 9D, (1.2) 


*HA has been supported in part by the NSF grant DMS-1521590. EO has been supported in 
part by CONICYT through FONDECYT project 3160201. AJS has been supported in part by NSF 
grant DMS-1418784. 

Department of Mathematical Sciences, George Mason University, Fairfax, VA 22030, USA. 
hantilOgmu.edu, 

^Departamento de Matematica, Universidad Tecnica Federico Santa Marfa, Valparaiso, Chile, 
enrique.otarolaSusm.cl 

^Department of Mathematics, University of Tennessee, Knoxville, TN 37996, USA. 
asalgadlOutk.edu 


1 



2 


H. Antil, E. Otarola, A.J. Salgado 


and the control constraints 


u € U^, (1.3) 

where is a nonempty, closed and convex subset of 12). The main source 

of difficulty and originality here is that the matrix A is not uniformly elliptic, but 
rather satisfies 


w(x)|eP <e-^(x)w(x)|Cp V^gR" (1.4) 

for almost every a: G 12. Since we allow the weight to vanish or blow up, this 
nonstandard ellipticity condition must be treated with the right functional setting. 
Problems such as arise when applying the so-called Caffarelli-Silvestre exten¬ 
sion for fractional diffusion when dealing with boundary control¬ 

lability of parabolic and hyperbolic degenerate equations [laiTiiisa and in the 
numerical approximation of elliptic problems involving measures [a ED]. In addi¬ 
tion, invoking Rubio de Francia’s extrapolation theorem [T9| Theorem 7.8] one can 
argue that this is a quite general PDF constrained optimization problem with an 
elliptic equation as state constraint, since there is no L^, only with weights. 

• Optimization -with point observations. Let Z C 12 with < oo. Given a 
set of prescribed values {yzjzGZ^ a regularization parameter A > 0, and the cost 
functional 


Jz{y,u) 




(1.5) 


the problem under consideration reads as follows: Find min Jz subject to 

— Ay = u in 12, y = 0 on 512, (1.6) 


and the control constraints 


u G U^:, (1.7) 

where Vz is a nonempty, closed and convex subset of L^(12). In contrast to standard 
elliptic PDF constrained optimization problems, the cost functional (II3|) involves 
point evaluations of the state. Note that these evaluations are not required for the 
state equation (fTDl) to be well posed. Additional assumptions must be made for 
this to make sense and, as will be seen below, the point evaluations of the state y 
in (II3|) lead to a subtle formulation of the adjoint problem. 

Problem (fFSt-lfTTl) finds relevance in numerous applications where the observa¬ 
tions are carried out at specific locations. For instance, in the so-called calibration 
problem with American options [I], in the optimal control of selective cooling of 
steel m, in the active control of sound [T] 148) and in the active control of vibrations 
ppl 155] ; see also [i [IDl EH El [S3j for other applications. 

The point observation terms in the cost (USD, tend to enforce the state y to have 
the fixed value j/z at the point z. Consequently, (IT31)~(fr7l) can be understood as 
a penalty version of a PDF constrained optimization problem where the state is 
constrained at a collection of points. We refer the reader to d Section 3.1] for a 
precise description of this connection and to [43] for the analysis and discretization 
of an optimal control problem with state constraints at a finite number of points. 


Weights, PDE constrained optimization, and error estimates 


3 


Despite its practical importance, to the best of our knowledge, there are only two 
references where the approximation of (IT31)~(IT:71) is addressed: [5] and m- In both 
works the key observation, and main source of difficulty, is that the adjoint state for 
this problem is only in with r G With this functional setting, 

the authors of [9] propose a fully discrete scheme which discretizes the control 
explicitly using piecewise linear elements. For n = 2, the authors obtain a 0{h) 
rate of convergence for the optimal control in the L^-norm provided the control and 
the state are discretized using meshes of size 0{h?) and 0{h), respectively; see [HI 
Theorem 5.1]. This condition immediately poses two challenges for implementation: 
First, it requires to keep track of the state and control on different meshes. Second, 
some sort of interpolation and projection between these meshes needs to be realized. 
In addition, the number of unknowns for the control is significantly higher, thus 
leading to a slow optimization solver. The authors of [5] were unable to extend these 
results to n = 3. Using the so-called variational discretization approach the 
control is implicitly discretized and the authors were able to prove that the control 
converges with rates 0{h) for n = 2 and for n = 3. In a similar fashion, 

the authors of m use the variational discretization concept to obtain an implicit 
discretization of the control and deduce rates of convergence of 0{h) and 0 (/i 2 ) 
for n = 2 and n = respectively. A residual-type a posteriori error estimator 
is introduced, and its reliability is proven. However, there is no analysis of the 
efficiency of the estimator. 

In Section |4] below we introduce a fully discrete scheme where we discretize the con¬ 
trol with piecewise constants; this leads to a smaller number of degrees of freedom 
for the control in comparison to the approach of [9]. We circumvent the difficul¬ 
ties associated with the adjoint state by working in a weighted H^-space and prove 
near optimal rates of convergence for the optimal control: 0{h\ log h\) for n = 2 and 
0{h'^ \ log/ip) for n = 3, respectively. In addition, we provide pointwise error esti¬ 
mates for the approximation of the state: 0{h\ log h\) for n = 2 and 0{h^\ log h\^) 
for n = 3. 

• Optimization with singular sources. Let V C ft he linearly ordered and with 
cardinality I = ^T> < oo. Given a desired state yd G and a regularization 

parameter A > 0, we define the cost functional 

Js{y, u) = \\\y- + ^||u||^,. (1.8) 

We shall be concerned with the following problem: Find min Js subject to 

— Ay = ^ UzSz in D, y = 0 on dft, (1.9) 

zei) 

where 6z is the Dirac delta at the point 2 and 

u = G U^, (1.10) 

where C K.* with U^, again, nonempty, closed and convex. Notice that since, 
for n > 1, 6z ^ H~^{fl), the solution y to (11.91) does not belong to H^{fl). Conse¬ 
quently, the analysis of the finite element method applied to such a problem is not 
standard [HI |50l EB] . We rely on the weighted Sobolev space setting described and 
analyzed in [50l Section 7.2]. 

The state (fTBl) . in a sense, is dual to the adjoint equation for dESD-dUD, where 
the adjoint equation has Dirac deltas on the right hand side. The optimization 


4 


H. Antil, E. Otarola, A.J. Salgado 


problem (EHD-dlll) is of relevance in applications where one can specify a control 
at a finitely many pre-specified points. For instance, references [71 |35] discuss 
applications within the context of the active control of sound and [HI |311 ISi] in the 
active control of vibrations; see also [24l [2^1^ . 

An analysis of problem (I1.8|) " (I1.10I) is presented in |^, where the authors use the 
variational discretization concept to derive error estimates. They show that the 
control converges with a rate of 0{h) and in two and three dimensions, 

respectively. Their technique is based on the fact that the state belongs to Wg^’'’(n) 
with r G {^ 2 ' addition, under the assumption that yd G L°°{Q) they 

improve their results and obtain, up to logarithmic factors, rates of 0{h?) and 0{K). 
Finally, we mention that [misT] study a PDF constrained optimization problem 
without control constraints, but where the controls is a regular Borel measure. 

In Section we present a fully discrete scheme for which we provide rates of con¬ 
vergence for the optimal control: in two dimensions and in three 

dimensions, where e > 0. We also present rates of convergence for the approxima¬ 
tion error in the state variable. 

Before we embark in further discussions, we must remark that while the intro¬ 
duction of a weight as a technical instrument does not seem to be completely new, 
the techniques that we use and the range of problems that we can tackle is. For 
instance, for integro-differential equations where the kernel g is weakly singular, the 
authors of m study the well-posedness of the problem in the weighted L^{g, (—r, 0)), 
space. Numerical approximations for this problem with the same functional setting 
were considered in [40] where convergence is shown but no rates are obtained. These 
ideas were extended to neutral delay-differential equations in [22l [23] where a weight 
is introduced in order to renorm the state space and obtain dissipativity of the un¬ 
derlying operator. In all these works, however, the weight is essentially assumed to 
be smooth and monotone except at the origin where it has an integrable singularity 
[miD] or at a finite number of points where it is allowed to have jump discontinu¬ 
ities 111123]. All these properties are used to obtain the aforementioned results. In 
contrast, our approach hinges only on the fact that the introduced weights belong 
to the Muckenhoupt class A 2 (see Definition 12. II below 1 and the pertinent facts from 
real and harmonic analysis and approximation theory that follow from this definition. 
Additionally we obtain convergence rates which, up to logarithmic factors, are opti¬ 
mal with respect to regularity. Finally we must point out that the class of problems 
we study is quite different from those considered in the references given above. 

Our presentation will be organized as follows. Notation and general considerations 
will be introduced in Section [2] Section [3| presents the analysis and discretization of 
problem Unii-iirsi). Problem (IT31)-(fr71) is studied in Section |T| The analysis of 
problem (II.81) -- (II.101) is presented in Section |5| Finally, in Section jO] we illustrate our 
theoretical developments with a series of numerical examples. 

2. Notation and preliminaries. Let us fix notation and the setting in which 
we will operate. In what follows D is a convex, open and bounded domain of K” 
(n > 1) with polytopal boundary. The handling of curved boundaries is somewhat 
standard, but leads to additional technicalities that will only obscure the main ideas 
we are trying to advance. By A < B we mean that there is a nonessential constant c 
such that A < cB. The value of this constant might change at each occurrence. 

2.1. Weights and weighted spaces. Throughout our discussion we call a 
weight a function w G such that oj{x) > 0 for a.e. x G M”. In particu- 
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lar we are interested in the so-called Muckenhoupt weights mm- 

Definition 2.1 (Mnckenhoupt class). Letr € (l,c)o) andw be a weight. We say 
that to G Ar if 

Cr,ui ■= sup a;(a:)dx^ dx^ < cx), 

where the supremum is taken over all balls B C K". 

From the fact that uj G Ar many fundamental consequences for analysis follow. 
For instance, the induced measure w dx is not only doubling, but also strong doubling 
(cf. [ini Proposition 2.2]). We introduce the weighted Lebesgue spaces 

n) = I?; e L°(fl) : J |x(x)|'’a;(x) dx < oo 

and note that [50l Proposition 2.3] shows that their elements are distributions, there¬ 
fore we can define weighted Sobolev spaces 

fl) = {i; e il) : G U{uj, ft) Vk : j^l < kj , 

which are complete, separable and smooth functions are dense in them (cf. |591 Propo¬ 
sition 2.1.2, Corollary 2.1.6]). We denote = lF^’^(a;,n). 

We define VL) as the closure of in fl) and set i?o (w, fl) = 

WQ^’^(a;, 11). On these spaces, the following Poincare inequality holds 

IklU^C^.n) < l|Vn|U.(^,n) Vn e (2.1) 

where the hidden constant is independent of v, depends on the diameter of O and 
depends on uj only through Cr,uj ■ 

The literature on the theory of Muckenhoupt weighted spaces is rather vast so we 
only refer the reader to [ninnii^ for further results. 

2.2. Finite element approximation of weighted spaces. Since, the spaces 
VF^’’’(a;,0) are separable for oj G Ar {r > 1), and smooth functions are dense, it is 
possible to develop a complete approximation theory using functions that are piecewise 
polynomial. This is essential, for instance, to analyze the numerical approximation 
of (11.21) with finite element techniques. Let us then recall the main results from [50] 
concerning this scenario. 

Let = {T} be a conforming triangulation (into simplices or n-rectangles) of 
n. We denote by T = {.^} a family of triangulations, which for simplicity we assume 
quasiuniform. The mesh size of £ T is denoted hy h^. Given G T we define 
the finite element space 

= {v^ G C°in) : e V{T), = 0} , (2.2) 

where, if T is a simplex, V{T) = Pi(T) — the space of polynomials of degree at 
most one. In the case that T is an n-rectangle V{T) = Qi(T) — the space of 
polynomials of degree at most one in each variable. Notice that, by construction, 
C IFg^’°“(n) C Wq’^{uj, n) for any r G (1, oo) and uj G Ar. 

The results of m show that there exists a quasi-interpolation operator 11^ : 
L^(fl) —>• which is based on local averages over stars and thus well defined for 
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functions in L^{n). This operator satisfies the following stability and approximation 
properties: 


Ik “ ^ ^5^lkllvEi.>-(w,n), 

\\^3rv\\w'^’^{u),n) ^ IkllvEC^(tj.n), 

Ik — ^ hs\\v\\w^’'-{ui,n), 


Vu e L^{oj,n), 
Vu e w^’''{uj,n), 
Vu G w^'^{uj,n), 
Vu G 


Finally, to approximate the PDF constrained optimization problems described in 
Section [T] we define the space of piecewise constants by 


U(^) = G : vs\T G Po(T)} . (2.3) 


2.3. Optimality conditions. To unify the analysis and discretization of the 
PDF constrained optimization problems introduced and motivated in Section [T] and 
thoroughly studied in subsequent sections, we introduce a general framework following 
the guidelines presented in [la [Ml Ea iMi 061 [58]. Let U and H be Hilbert spaces 
denoting the so-called control and observation spaces, respectively. We introduce the 
state trial and test spaces Yi and Xi, and the corresponding adjoint trial and test 
spaces Y 2 and X 2 , which we assume to be Hilbert. In addition, we introduce: 

(a) A bilinear form a : (Yi + Y 2 ) x (Xi -|- X 2 ) ^ M which, when restricted to either 
Yi X Xi or Y 2 X X 2 , satisfies the conditions of the BNB theorem; see [Ml Theorem 
2 . 6 ]. 

(b) A bilinear form 6 : U x (Xi +X 2 ) —>■ M which, when restricted to either U x Xi or 
U X X 2 is bounded. The bilinear forms a and b will be used to describe the state 
and adjoint equations. 

(c) An observation map C : Dom(C') C Yi -|- Y 2 ^ H, which we assume linear. 

(d) A desired state yd G H. 

(e) A regularization parameter A > 0 and a cost functional 

Y,xV5{y,u)^Jiy,u) = ^\\Cy-yd\g+^\\u\\l (2.4) 

All our problems of interest can be cast as follows. Find min J subject to: 

?/G Yi : a{y,v) = b{u,v) Vu G Xi, (2-5) 


and the constraints 


U G Uad; 


( 2 . 6 ) 


where Uad C U is nonempty, bounded, closed and convex. We introduce the control 
to state map S : U Yi which to a given control, u G U, associates a unique state, 
y{u) = Su G Yi, that solves the state equation ()2.5I) . As a consequence of (jlj) and 
©, the map S' is a bounded and linear operator. With this operator at hand we can 
eliminate the state variable y from (12.41) and introduce the reduced cost functional 

U 9 u !-)■ j{u) = -IICSu — ydWm + '^ll^llu- (2-7) 

Then, our problem can be cast as: Find minj over Uad- As described in (gj) we have 
A > 0 so that j is strictly convex. In addition, Uad is weakly sequentially compact in 
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U. Consequently, standard arguments yield existence and uniqueness of a minimizer 
[551 Theorem 2.14]. In addition, the optimal control u G Uad can be characterized by 
the variational inequality 


j'(u)(M-u)>0 VuGUad, 

where j'{w) denotes the Gateaux derivative of j at w [58l Lemma 2.21]. This varia¬ 
tional inequality can be equivalently written as 

6 (u — u, p) -I- A(u, u — u)u > 0 Vu G Uad, (2.8) 

where p = p(u) denotes the optimal adjoint state and solves 

p G X 2 : a{v,p) = {Cy-yd,Cv)u Vu G ¥ 2 . (2.9) 

The optimal state y = y(u) G Yi is the solution to (12.51) with u = u. 

2.4. Discretization of PDE constrained optimization problems. Let us 

now, in the abstract setting of Section [2^ study the discretization of problem (12.41) - 
dSll). Since our ultimate objective is to approximate the problems described in Section 
|T] with finite element methods, we will study the discretization of (I2.4l) " (l2.6p with 
Galerkin-like techniques. 

Let /i > 0 be a parameter and assume that, for every h > 0, we have at hand 
finite dimensional spaces C U, C Xi, X2 C X2, ¥]* C Yi and Y2 C Y2. We 
define U^^ = U^ D Uad, which we assume nonempty. About the pairs (X]*, Y(*), for 
i = 1,2, we assume that they are such that a satisfies a BNB condition uniformly in 
h; see [501 §2-2.3]. In this setting, the discrete counterpart of (I2.4l) - (l2.6p reads: Find 

min J{yh,Uh) ( 2 - 10 ) 

subject to the discrete state equation 

y/iGY^: a{yh,Vh) = b{uh,Vh) G X^, (2.11) 

and the discrete constraints 


un G U^d- (2.12) 

As in the continuous case, we introduce the discrete control to state operator Sh, 
which to a discrete control, Uh G U;,, associates a unique discrete state, yh = yh(uh) = 
ShUh, that solves ( 12 .lip . Sh is a bounded and linear operator. 

The pair {yh,^h) G Y(* x U^^j is optimal for (I2.10I) - (I2.12I) if y^ = y/,(u/j) solves 
(I2.1ip and the discrete control D^, satisfies the variational inequality 

jhi^h)iuh-Uh) >0 Vu,,GU^d, 


or, equivalently, 


b{uh - Uh,Ph) + - Uh)v >0 Vit/, G Uad, (2.13) 

where the discrete adjoint variable pu = P/t(uft) solves 

Ph€X2' a{vh,Ph) = (Cyh - yd,Cvh)M G ¥ 3 . (2.14) 
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To develop an error analysis for the discrete problem described above, we intro- 
duce Hu the U-orthogonal projection onto U^. We assume that IluUad C U^^j. In 
addition, we introduce two auxiliary states that will play an important role in the 
discussion that follows. We define 


a{yh,Vh) = b{u,Vh) Vu/i e Xj", (2.15) 

i.e., Yh is defined as the solution to (12.111) with Uh replaced by u. We also define 

p?, G X 2 : a{vh, ^h) = (Cyh - 2/d, Cvh)M '^Vh G Y 2 , (2.16) 

this is, ph is the solution to (12.1411 with yh replaced by yn- 

The main error estimate with this level of abstraction reads as follows. 

Lemma 2.2 (abstract error estimate). Let (y, u) G Yi x Uad and {yh,Uh) G 
X be the continuous and discrete optimal pairs that solve (12.41) - (12.61) and 
(12.101) - (12.121) . respectively. If 

GXjnx^, y^-y^eYjnY 2 \ ( 2 . 17 ) 


then, for any e > 0 , we have the estimate 

||u - < Cl (Up - Pdlll, +j'(u)(nuu - u) + IlHuu - u||^) 

/ _ . \ 2 
{C{yh — Yh), Cvp)m \ 


+ e 


sup 

i>peY& 


Ikplks 


(2.18) 


where the constant ci depends on and but does not depend on h. 

Proof Since by definition Uad U Uad and by assumption IluUad C U^d, we set 
u = Ufi and Uh = Huu in (12.81) and (12.131) respectively. Adding the ensuing inequalities 
we obtain 


A||u - u,,||^ < b{uh - u, p - p,,) + 6(nuu - u, p,,) + A(u/j,nuu - u)u. (2.19) 

Denote I = b{Uh — u, p — p;,). In order to estimate this term we add and subtract 
Ph to obtain 


I = b{uh - u, p - pi,) + b{uh -w,Ph- Pd). (2.20) 

Since pd is the unique solution to (12.161) . we have that 

a(z;p, pd - Pd) = ((^(vd - yd),C'up)H Vwp G Y^. (2.21) 

Similarly, since yd solves (12.151) . we derive 

a(yd - yd,Uy) = 6(ud - u,Uy) Vuy G Xk 

Set Vp = Yh — Yh and Vy = ph — Ph which, by assumption (12.1711 . are admissible test 
functions to obtain 

6 (nd -U,ph- Ph) = {C{yh - yd), C{% - yd))H < 0. 

This, and the continuity of the bilinear form b allow us to bound (12.201) as follows: 

I < b{Uh - u,p- p,,) < ^||u - u,.||^+ li^llp- p,,|||^. 
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where ||&|| denotes the norm of the bilinear form b. 

Let us now analyze the remaining terms in (j2.19|) . which we denote by II. To do 
this, we rewrite II as follows: 


II = 6(nuu - u, p) + A(u, IIuu - u)i[j + A(u,i - u, IIuu - u)u 

+ 6(nuu - u, - p) + 6(nuu -u,ph- ph). 

Now, notice that 

6(nuu - u, p) + A(u, IIuu - u)u = /(u)(nuu - u) 

and 

x{uh - u,niuu - u)u < ^||a - + i||o - npuHu. 

Next, since the bilinear form b is continuous, we arrive at 


6(nuu - u,p/, - p) < ■^||nuu-u||^ + -^||p-p^|||^. 

The remaining term, which we will denote by III, is treated by using, again, that the 
bilinear form b is continuous. This implies that for any e > 0 


III := &(nuu -u,ph- M < -^IlHuu - u||f, + ||6|||||p/. - p/.|||^. 

From (12.211) and the fact that the discrete spaces satisfy a discrete BNB condition 
uniformly in h we conclude 


\Ph - Ph\ 


< 


sup 


jCjyh — 9h), CVp)m 

Ikplkj 


Collecting these derived estimates we bound the term II. 

The estimates we obtained for I, II and III readily yield the claimed result. □ 
The use of this simple result will be illustrated in the following sections. 
Remark 2.3 (on the choice of e). In (I2.18p the choice of e can be arbitrary. A 
judicious choice will be fundamental in the error analysis of the optimization problem 
with point observations (lEsii-dni). 

3. Optimization with nonuniformly elliptic equations. In this section we 
study the problem (in])-(IE3l) under the abstract framework developed in Section 
12.31 Let n C M" be a convex polytope («■>!) and w € A 2 (R") where the A 2 - 
Muchenkhoupt class is given by Definition 12.II In addition, we assume that A : D —>■ 
M" is symmetric and satisfies the nonuniform ellipticity condition (11.41) . 

3.1. Analysis. Owing to the fact that the diffusion matrix A satisfies (|1.4|1 with 
u) £ A 2 (K”), as shown in [ 21 ] . the state equation ( 1121 ) is well posed in iLg(a;,D), 
whenever u £ , D). For this reason, we set: 

• m = L^{uj,n) and (7 = id. 

• u = 

• Xi = X 2 = Yi = ¥2 = Hq{uj, D), and 

a{vi,V2) = / yv2{xyA{x)Vvi{x)dx, 

Jo. 

which, as a consequence of (11.41) with w £ A 2 (M") and the Poincare inequality (12.11) . 
is bounded, symmetric and coercive in iLg (w, D). 
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• Notice that, if Vi G L'^{w ^,ri) and V 2 G then 

b{vi,V2) = {vi,V2)L^(n) < ||fi||L2(w-i.n)lk2||L2(,^_n) ^ lki||L2(,^-i_n)||V?;2||L2(,^_n), 

where we have used the Poincare inequality (EH). 

• The cost functional as in EH- 

For a, b G M, a < b we define the set of admissible controls by 

= {w G 12) : a < m < b a.e. x G 12} , (3.1) 

which is closed, bounded and convex in L^(a;“^,12). In addition, since A > 0 the 
functional (EH is strictly convex. Consequently, the optimization problem with 
nonuniformly elliptic state equation (EH-(EH has a unique optimal pair (y, u) G 
_ffQ(w,12) X L^(w“^,12) [58l Theorem 2.14]. In this setting, the first order necessary 
and sufficient optimality condition (12.81) reads 


(P,M-u)i 2 (o) + A(u,'u-u)i 2 (^-i fj) > 0 VmGU^, (3.2) 

where the optimal state y = y(u) G Hq{u!, 12) solves 

a(y, v) = (0,w)L 2 (n) Vu G Hg{uj, 12) (3.3) 

and the optimal adjoint state p = p(u) G Hq{uj, 12) solves 

a{v,p) = {y-yd,v)L2(u:,n) V?; G i?}(a;, 12). (3.4) 

The results of El]: again, yield that the adjoint problem is well posed. 

3.2. Discretization. Let us now propose a discretization for problem (EH- 
EH, and derive a priori error estimates based on the results of Section 12.41 Given a 
family T = {,^} of quasi-uniform triangulations of 12 we set: 

• = 1[J(1F), where the discrete space 'U(^) is defined in (12.31) . 

• ^ad = n U.A, where the set of admissible controls is defined in (13.11) . 

• Ho is the L^(a;“^, 12)-orthogonal projection onto which we denote by n^^-i 

and is defined by 

(n,^-i'^)|r = _ ^ [ uj~^ {x)v{x) dx (3.5) 

^ [T) 

The definition of yields that II^^-iU^ C 

• = X 2 = = Y 2 = V(.^), where the discrete space V(fY) is defined in (EH- 

Notice that, since X(* = X 2 = Y^ = Y 2 , assumption (I2.17|) is trivially satisfied. 

We obtain the following a priori error estimate. 

Corollary 3.1 (a priori error estimate). Let 0 and Uu be the continuous and 
discrete optimal controls, respectively, //y, p G then 

||U - U/,||i2(,^-qo) < ||U -n^-iu|E2(,^-qn) + II^P “ (^p) ||L2(^-qn) 

+ ^5"(l|y||ff2(^_n) + l|p||_f/2(,^_o)), 

where the hidden constant is independent of h^r- 

Proof. We invoke Lemma [2.21 and bound each one of the terms in (I2.18p . First, 
since y, p G 12), the results of [50] imply that 

Up — P/t||ffi(a;,n) ^ hsr (||y||ff2(a,,n) + I|p||ff2(;^_f2)) ■ 
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Indeed, since p solves (13.41) and Ph solves (I2.16P , the term p — Ph satisfies 
a{vh, P-Ph) = {y- 9h, Vh)L^u;.n) ^Vh e 

Adding and subtracting the terms and p appropriately, where 11 5 ^ denotes the 
interpolation operator described in H2.21 and using the coercivity of a we arrive at 

Up - Ph\\Hi(u;,Q) < Up - ^ 3 rp\\Hi(u;,Q) + lly - 

Using the regularity of p and y we obtain the claimed bound. 

We now handle the second term involving the derivative of the reduced cost j. 
Since it can be equivalently written using (12.81) . invoking the definition of n^^-i given 
by (13.Sp . we obtain 

/(u)(n^-iu - D) = (p,n^-iu - D)L 2 (n) + A( 0 ,n^-i 0 - 0 )i 2 (^-yn) 

= (wp- n^-i(a;p),n^-iu - D)i2(^-yo) - A||n^-in - 

< ||wp — n^^-i (wp)||^ 2 (^-yn) + 

The Poincare inequality (12.11) . in conjunction with the stability of the discrete 
state equation ( 12 . 111 ) . yield 

iVh -yh,Vh)L'^(u;,a) < ||y?i - y?i||ffi(c^,n)||w/i||ffi(a;.n) 

^ l|u - u/t||L 2 (^-yn)lk?i||ffi(a;,n)i 

for all Vh €Y(3^). This yields control of the last term in (12.181) . 

These bounds yield the result. □ 

Remark 3.2 (regularity of y and p). The results of Corollarv l3.1l relv on the fact 
that y, p G Reference [16] provides sufficient conditions for this to hold. 

Theorem 3.3 (rate of convergence). In the setting of Corollo,rv \3.1\ if we addi¬ 
tionally assume that wp € then, we have the optimal error estimate 

l|u - u/i||L2(^-yn) < h^r (||y||H2(cj,n) + ||p||H2(i.d.n) + ll‘^p||Hi(i.j-i.n) + llullHqw-yn)) ) 

where the hidden constant is independent of h^. 

Proof We bound ||u - n^-iu||^ 2 (^-yn) and ||wp - n^-i(wp)||i 2 (^-qn)- Using 
that ujp € , n) and a Poincare-type inequality (SO] Theorem 6.2], we obtain 

l|wp — n(^-i(wp)||L2(,^-qQ) < /l5"||wp||//i(^-qn)- 

Now, to estimate the term IIj^-iD — u, it is essential to understand the regularity 
properties of 0 . From [SS] Section 3.6.3], u solves (13.2p if and only if 

u = max a, min 

The assumption cup G iy^(cu“^,n) immediately yields 0 G ,Vl) Theorem 

A.l], which allows us to derive the estimate 





Collecting the derived results we arrive at the desired estimate. □ 
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4. Optimization with point observations. Here we consider problem 
(ini). Let n C M" be a convex polytope, with n S {2,3}. We recall that Z C H 
denotes the set of observable points with < oo. 

4.1. Analysis. To analyze problem (ll.5l) " (ll.7l) using the framework of weighted 
spaces we must begin by defining a suitable weight. Since < oo, we know that 
dz = min{|z — z'\ : z,z' & Z, z ^ z'j > 0. For each z € Z we then define 


and the weight 


dz(a;) 



Wz{x) 


log"d,(x) 


{ zuz{x), 3z e Z : d^{x) < -, 

22 —n 

- TT -, otherwise. 

log' 2 ’ 


(4.1) 


As [Sni Lemma 7.5] shows, with this definition we have that ccj S A 2 . With this 
H 2 -weight at hand we set: 

• H = and C = e zSzi where {ez\zez is the canonical basis of H. 

• U = L2(rj). 

. Xi = Yi = 

• X 2 = HqIw, n) and Y 2 = Hq{zu~^, n) and 


a{v,w) = / Vv{xy ■ \/w{x) dx, 
Jn 


which is bounded, symmetric and coercive in and satisfies the conditions of 

the BNB theorem in iLg (ccj, H) x iLg ( 07 “^, ft) [SOI Lemma 7.7]. 

• 6 (-, •) = Tlio results of [SUl Lemma 7.6] guarantee that, for n < 4, the 

embedding fl) ^ holds. Therefore, 


byuV2) < lki||L2(n)||'C2||ffi(z^.n)- 

For a, b G R with a < b we define the set of admissible controls by 

Vz = {w G : a < u < b, a.e. a; G 12} . (4.2) 

With this notation, the pair (y, u) G Hg (^) ^ ^^(17) is optimal for problem (11.51) - 
(ini) if and only if y solves 

y G Hg{n) : a{y,w) = {u,w)L 2 {n) Vw G i7{(17), (4.3) 

and the optimal control u satisfies 

(p,M - a)i, 2 (n) + A(Q,u-u)L 2 (n) > 0 Vm G U^, (4.4) 

where the adjoint variable p G Hyw, 17) satisfies, for every w G i7}(z77“^, 17), 

(4.5) 


Notice that, since 17 is a convex polytope and n < 4, we have that y G 77^(17) ^ ^(17), 
so point evaluations are meaningful. Therefore, since Sz G i7{('cc7,17)' for z G 17, using 
[SOI Lemma 7.7] we have that the adjoint problem is well posed. 
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4.2. Discretization. For a family T = {^} of quasi-uniform meshes of we 
set: 

# = U(^), where is defined in (12.31) and nU^, where Vz is 

defined in (14.21) . The operator IIu = 11^2 is the standard L^(n)-projection: 

(n^ 2 u)| 7 ’ = J' v{x) da; VT S 

• Xf = 

To shorten the exposition we define 

o-jr =log/i.y|. (4.6) 

With this notation, the error estimate for the approximation (I2.10I) - (I2.12I) to 
problem (fra-ifTTi) reads as follows. 

Corollary 4.1 (a priori error estimates). Let u and Uh be the continuous and 
discrete optimal controls, respectively. Assume that is sufficiently small. Ifn = 2, 
then we have 

||u — u/j||i 2 (Q) < ||u — ni 2 u|| 2 , 2 (Q) -I- a^r (||Vp||L 2 (CT,n) + II Vy||Loo(Q)) . (4.7) 

If n = 3, we have 

||u - u,,||i 2 (Q) <|log/i.y|||u -ni 2 u||i 2 (Q) +as (||Vp||l 2 (^_o) -f ||Vy||Loo(n)) , (4.8) 

where is defined in (gH) and the hidden constants are independent of SA. 

Proof. We follow Lemma 12.21 with slight modifications. The term I in (12.201) is 
estimated in two steps. In fact, since (u^, — u,ph — Ph)L^{Q) < 0, we have 

I < (u?, - U, p - Mlyq) < ^I|u/i - 0||i2(o) + ^IIp - PhWh^n)- 
We now analyze the second term of the previous expression. Let G Y{.tY) solve 

a{wh,qh) = '^iy{z) - yz)wh{z) 'iwhGY{^). (4.9) 

The conclusion of [50l Corollary 7.9] immediately yields 

Up - qhllL2(o) < crjrl|Vp||i2(^^n), 

so that it remains to estimate qh — ph- We now invoke [501 Theorem 6.1] with p = q = 
2, p = 1 and ui = w. Under this setting the compatibility condition [50l inequality 
(6.2)] is satisfied, and then (SHl Theorem 6.1] yields 

||q?i - P/i||L 2 (n) < ||V(q/i - P/i)||l 2 (^^o)_ 

where the hidden constant depends on U and the quotient between the radii of the 
balls inscribed and circumscribed in U. Since q/j solves (14.91) . the discrete inf-sup 
conditions of (Si Lemma 7.8] and the fact that bz G IIo(zi7 yield 


l|v(q,, - < l|y-y/«(u)||L=o(o). 
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We now recall that % is the Galerkin projection of y. In addition, since n G {2,3}, 
n is a convex polytope and u G we have that y G (cf. [^1^ ITT] !. 

Therefore standard pointwise estimates for finite elements [551 Theorem 5.1] yield 

l|y-y?i||L“(n) < ^.y|logh^|||Vy||i=o(n). (4.10) 


In conclusion. 


I < “ u||^2(Q) + calr (||Vp||^2(^_f2) + ll^yllL“(n)) > 

for some nonessential constant c. 

We estimate the term /( 0 )(ni 2 n — u) as follows: 

/(u)(nL2n - u) = (p + Au,ni2u - a)L2(n) = (p + Au - ni2(p + AD),ni2u - u)i2(Q) 

< ~ ^Ili 2 (n) + 2 Up “ n 2 , 2 p||| 2 (Q) 

< 2 ~ ^Ili2(n) + ccr|r||Vp|||2(^_f2)) 

for some nonessential constant c. We have used the properties of 11^2, together with 
the Sobolev-Poincare inequality of [501 Theorem 6.2]; see also [SOI Corollary 7.9]. 

We now proceed to estimate the term III in the proof of Lemma [2.21 as follows: 

III := 5(11^25 - u, p/j - p/i(u)) = (11^25 -u,ph- Ph{u))L^(Q) 

1 s 

< —||ni2U - + -\\ph - Ph(u) - ^L^Ph - P/i(u))||i2(n) 

- “ °lli=(n) + cecr^llV(p,, - P;i(u))lli2(^_n), 

where we have used the properties of 11^2 together with the Sobolev-Poincare in¬ 
equality of m Theorem 6.2]. In the previous estimate £ > 0 is arbitrary and c is 
a nonessential constant; see Lemma [2.21 Using now the fact that Sz G fl)', 

and the discrete inf-sup stability of [50l Lemma 7.8], we have that 

l|P/i - P?i||ffi(ro.n) ^ lly?i - y?i||L“(n) 

< ^sWyh - y/i|lHi(n) < i^^lju - u/tllL2(o), (4.11) 

where i^^ is the mesh-dependent factor in the inverse inequality between L°°(n) and 
H^{Q) (see [S] Lemma 4.9.2] for n = 2 and [20j Lemma 1.142] for n = 3): 

= (1 -I- 1 log h^-j) 2 if n = 2, and i^r = if n = 3. (4.12) 

Therefore, we have derived 

III ^ ~ ^111.2(0) + cecr|ri|r||u — U/ijj^2(Q). 

We now examine the product for sufficiently small: 

if n = 2, 

\\\oghs\ if n = 3. 
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Collecting the derived estimates we arrive at the desired estimates (I 121 )-gll) by 
considering a judicious choice of e. To be precise, when n = 2, we set e = 1. If n = 3, 
we consider e = c| log hsr\~'^- 0 

Proposition 4.2 (regularity of u). 7/0 solves (11.51) - (11.71) then u G 
Proof. From [58l Section 3.6.3], 0 solves (|4.4I) if and only if 


u = max 


■| a, min 



This immediately yields u G 77^(ro,n) by invoking Theorem A.lj. □ 

Using this smoothness and an interpolation theorem between weighted spaces 
we can bound the projection error in Corollary 14.11 and finish the error estimates 
dllD-gS as follows. 

Theorem 4.3 (rates of convergence). In the setting of Corollary we have 

l|u - a?,||i 2 (n) < h3r\\ogh^\ (||Vp||i 2 (^_n) + ||Vy||i.»(n)) (4.13) 

for n = 2. Ifn = 3, we have 

||u - u,,||i 2 (o) < h^^\ log/ij7p (||Vp||i 2 (^_n) + l|Vy||Loo(o)) . (4.14) 

The hidden constants in both estimates are independent of , the continuous and 
discrete optimal pairs. 

Proof. We only need to bound the projection error jju — 11^2 011 ^ 2 ( 0 ). Proposi¬ 
tion yields u G then, invoking [SHI Theorem 6.2], we derive 

||u — 11 ^ 2011 ^ 2 ( 0 ) < cr 5 "||Vu]ji 2 (j,,_o)- 

Substituting the previous estimate in the conclusion of Corollary 14.11 we derive the 
claimed convergence rates. □ 

On the basis of the previous results we now derive an error estimate for the 
approximation of the state variable. 

Theorem 4.4 (rates of convergence). In the setting of Corollary we have 

l|y - yh||L“(n) < hs\\oghs\ (||Vp| 1 l 2 (^_o) + llVyl]Loo(o)) , (4.15) 

for n = 2. If n = 3, we have 

l|y - yh||L“>(n) < log/ij 7 p (||Vplli, 2 (^_o) + ||VyllL=o( 0 )) . (4.16) 

The hidden constants in both estimates are independent of ST, and the continuous and 
discrete optimal pairs. 

Proof. We start with a simple application of the triangle inequality: 


jjy- yh||Lo°(n) < lly -ylU-cn) + lly -yh||Lo°(n), (4.17) 

where y solves a(y, v) = (0^, v) for all v G 77(J(0). The second term on the right hand 
side of the previous inequality is controlled in view of standard pointwise estimates 
for finite elements. In fact, [5S1 Theorem 5.1] yields 


||y-y/i||L~(n) < ^j7|log/i.ylllVyllLcx,(o). 


(4.18) 
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To control the first term on the right hand side of (I4.17P , we invoke the results of [TT] 
and conclude that there is r > n such that 

l|y - y||iyi'r(n) < ||u - 

This, in conjunction with the embedding 0(^1) for r > n, yields 

l|y-y||L“(n) ^ l|u - unh^Q)- 

In view of (|4.18|) , the previous estimate and the results of Theorem 14.31 allow us to 
derive the desired error estimates. □ 

5. Optimization with singular sources. Let us remark that, since the for¬ 
mulation of the adjoint problem (|4.5I1 led to an elliptic problem with Dirac deltas on 
the right hand side, the problem with point sources on the state (ll.Sp - dl.lOl) is, in a 
sense, dual to one with point observations (fra-ifTTi). In the latter, the functional 
space for the adjoint variable is the one needed for the state variable in (I 1 . 8 l) - p.l 0 l) . 
The analysis will follow the one presented in Section [4.2l It is important to comment 
that problem (I1.8l) - (ll.l0p has been analyzed before. We refer the reader to |30] for 
the elliptic case and to [28l Ell EH [57] for the parabolic one. It is our desire in this 
section to show how the theory of Muckenhoupt weights can be used to analyze and 
approximate problem (11.81) ^ (11.101) . In doing this, it will be essential to assume that 
dist((D, dil) > dx> > 0. Set 

• H = L^{n) and C = id. 

• U = K'. 

• Yi = Hq{zu, fl) and Xi = Hq{zu ~^with w defined, as in Section ElTl by (14.11) . 

• Y 2 = X 2 = and 

a{v,w) = / Vn(a;)^ Vr(;(a;) dec. 

Jq 

• The bilinear form 6 : U x (Xi -I- X 2 ) is 

b{-v,w) = 

zeT) 

Since, for z € fl, Sz € fl)', we have that b is continuous on M} xHq{w~^ , fl). 

For a, b G R* with we define the set of admissible controls as 

Ui = {u G : az < Mz < bz, Vz G V} . 

The space of controls is already discrete, so we set = U and = U 5 . Finally 
we set, for i = 1, 2, X^ = Y^ = Y{3^). Since the bilinear form b is not continuous on 
U X X 2 , we need to slightly modify the arguments of Lemma EHl In what follows, for 
V G (7(11) and w G R^ we define 

{v,w)t) '^v{z)wz- (5.1) 

zeT> 

In this setting, the main error estimate for problem (jl. 8 l) - (ll.lOI) is provided be¬ 
low. We comment that our proof is inspired in the arguments developed in m 
Theorem 3.7] and [52] . 

Theorem 5.1 (rates of convergence). Letu anduh be the eontinuous and discrete 
optimal eontrols, respectively, and assume that for every q G (2, 00 ), yd G Let 

e > 0 and fli be such that 22 (s Di (s D. Ifn = 2, then 

l|u - U^iIIrI < (||u||r1 -I- ||p||h2(Q) -I- ||p||w2,r(ni)) . 


(5.2) 
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On the other hand, if n = S, then 

l|u - Uh\\Ri < (||u||ri + ||p||ff2(n) + ||p||iv2,r(ni)) , (5.3) 

where r < n/(n — 2). The hidden constants in both estimates are independent of ^, 
and the continuous and discrete optimal pairs. 

Proof. We start the proof by noticing that, since y — yd & Lf{Vl) and is convex, 
standard regularity arguments m yield p S ^ C{fl) . This guarantees that 

pointwise evaluations of p are well defined. Moreover, since, in this setting, = Uad 
estimate (j2.19ll reduces to 

A||u - u?,||gi < (p - Ph, - u)x., 

where (•, •)-d is defined in (15.ip . Adding and subtracting the solution to (12.161) ph, we 
obtain that 


A||u - < (p - Ph, Uh - u)t) + {ph - Ph, Uh - u)-D. (5.4) 

This, in view of {ph - Ph, W - u)x> = -||y?i - y/i|li 2 (Q), implies that 

M\^ - wIIr; + \\yh - yhWh^Q) < (p - h, w - u)x) 

< (p - Ph, Uh - u)x> + (q/t - Ph, W - u)x>, (5.5) 

where is defined as the unique solution to 

a{wh,Ph) = {y - yd,Wh)mn) \lwh&y{^). 

Since, by assumption, we have that d-p > 0 we can conclude that there are 
smooth subdomains flo and fli such that (D C flo ^ In view of dESl), this 

key property will allow us to derive interior L°“-estimates for p — q/t and q^, — ph- 

Let us first bound the difference p — q/j. To do this, we notice that, since y G 
WQ^’'*(n) for s < n/{n — 1), a standard Sobolev embedding result implies that y G 
L’"(n) with r < ns/{n — s) < n/(n —2). Then, on the basis of the fact that yd G L^{Tl) 
for q < oo, interior regularity results guarantee that p G W^’'’(r2i) for r < nl{n — 2). 
Consequently, since q/j corresponds to the Galerkin approximation of p, [541 Theorem 
5.1] yields, when n = 2, that for any e > 0 we have 

Up - q?i||L“(ao) ^ (^^''Ilp|lw2>'-(ni) + ^^l|p||H2(a)) . (5.6) 

When n = 3, we have that p G Hq{Q) (~l for r < 3 and, as a consequence. 

Up - qft||L“(no) ^ ''IIPllw2>'-(ni) + ^^||p||H2(f2)) . (5.7) 

It remains then to estimate the difference Ph = Ph — Ph- To do so we employ a 
duality argument that combines the ideas of [50j Corollary 7.9] and [44l|52]. We start 
by defining ip G Tfg(Q) as the solution to 

a{v,ip)= [ sgn(y-y?,)u Vu G i7o(0), (5.8) 

Jn 

where % solves (12.151) . Notice that Jj sgn(y — y/i)jjicx=(Q) < 1 for all gT. Therefore, 
[Ml Theorem 5.1] followed by [S^ Theorem 5.1] lead to (see also [301 Lemma 3.2]) 

\\p - Ph\\L--(Qo) ^hl^l'^Oghsl^ 


(5.9) 
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where Lph is the Galerkin projection of Lp and the hidden constant does not depend 
on ^ nor Lp. In addition we have that (p G n Hq{^) ^ Hq{w~^^VL) [5(7[ 

Lemma 7.6]. Therefore (/? is a valid test function in the variational problem that y 
solves. Then, using the continuity of the bilinear form a and Galerkin orthogonality 
we arrive at 

l|y-y/«llLi(f2) = / sgn{y-yh){y-9h) = a{y 

Jn 

= a{y, p-<ph) = {p- Ph, u)-D < ||u||Ri \\p - 7’/i||L“(no) 

< h|r|log/lj7p||u||Ri, 

where in the last step we used estimate (15.91) . 

We now recall that Pu solves 

a{wh., Ph) = (y - Wh)L^n) ' 9 wh e ^i^), 

an inverse inequality and a stability estimate for the problem above yield 

ll-P'»llL“(n) ^ < i^lly “ y?i||Li(n) ll^?i||L“(n)i 

where is defined in (14.121) . In conclusion, 

ll^’/i||L“(n) ^ ^s\\y -9h\\mn) < i|^^^|loghj 7 p||u||Ri. (5.10) 

Combining the obtained pointwise bounds for p — q/j and — ph we obtain the 
desired estimates. □ 

Remark 5.2 (comparison with the literature). Reference [30] claims to obtain 
better rates than those in Theorem 15.11 namely they can trade the term by a 
logarithmic factor \\ogh^\^ with s > 1 but small. However, when following the 
arguments that lead to this estimate (see [SO] formula (3.40)]) one realizes that a 
slight inaccuracy takes place. Namely, the authors claim that, for s <n/(ri— 1), 

I log h1 < I log h1. 

However, 3 — n/s < 4 — n which, for n = 2 or n = 3 reduces to the estimates that we 
obtained in Theorem 15. II 

To conclude, we present an error estimate for the state variable. 

Corollary 5.3 (rates of convergence). In the setting of Theorem \5. 1\ we have 

l|y - yh\\L^(n) < (IjujlRi + ||p||ff2(o) + llpllvE2.>-(ni) + ||Vy|lL2(^_f2)) . 

The hidden constant is independent of ST and the continuous and discrete optimal 
pairs. 

Proof. A simple application of the triangle inequality yields 

l|y - y/i||L2(n) < jjy - yh\\L^{n) + ||y?i - yh\\L^{n), ( 5 -ii) 

where yt solves (12.151) . To estimate the first term on the right-hand-side of the 
previous expression, we invoke [SOI Corollary 7.9] and arrive at 

l|y - y/i||L2(n) < crj 7 |lVy]lL 2 (^ O). 

Using (|531) and the results of Theorem Q] we bound the second term on the right 
hand side of (|5.1ip . This concludes the proof. □ 
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6. Numerical Experiments. In this section we conduct a series of numerical 
experiments that illustrate the performance of scheme (I2.10I) - (I2.12I) when is used to 
approximate the solution to the optimization problem with point observations studied 
in Section |4] and the one with singular sources analyzed in Section [5l 

6.1. Implementation. All the numerical experiments that will be presented 
have been carried out with the help of a code that is implemented using C++. The 
matrices involved in the computations have been assembled exactly, while the right 
hand sides and the approximation error are computed by a quadrature formula which 
is exact for polynomials of degree 19 for two dimensional domains and degree 14 
for three dimensional domains. The corresponding linear systems are solved using 
the multifrontal massively parallel sparse direct solver (MUMPS) [31 0]. To solve 
the minimization problem (I2.10p - (I2.12I) we use a Newton-type primal-dual active set 
strategy |S31 Section 2.12.4]. 

For all our numerical examples we consider A = 1 and construct exact solutions 
based on the fundamental solutions for the Laplace operator: 


cj){x) 


zee 

— y ^ , 

47r ^ \x — zV 
zee ' ' 


if U C 

if U C 


( 6 . 1 ) 


where, depending on the problem, & = Z ov @ = 7). 

We must remark that the introduction of weights is only to simplify the analysis 
and that these are never used in the implementation. This greatly simplifies it and 
allows for the use of existing codes. 

6.2. Optimization with point observations: one point. We set n = 2 and 

U = (0,1)^. We consider the control bounds that define the set U^: as a = —0.4 and 
b = —0.2. To construct an exact solution to the optimization problem with point 
observations, we slightly modify the corresponding state equation by adding a forcing 
term f G L^(Sl), i.e., we replace (|4.3I) by the following problem: 

yGLfo(U): o(y, in) = (u-|-f, w) Vin G iLo(U). (6.2) 


We then define the exact optimal state, the observation set and the desired point 
value as follows: 


y(a;i,X2) = 3211x2(1 - a;i)(l - 0:2), Z = {( 0 . 5 , 0 . 5 )}, y(o.5,o.5) = 1 - 

The exact optimal adjoint state is given by (EH) and the right hand side f is computed 
accordingly. 

To present the performance of the fully discrete scheme ()2.10|) - ()2.12L we consider 
a family of quasi-uniform meshes {^k}k=i- We set N{k) = that is, the total 

number of degrees of freedom. In addition, we denote by EOCq(fc) the corresponding 
experimental order of convergence associated to the variable q, which is computed 
using 


EOC{k) = 


ln(eq(fc - l)/eq(fc)) 

\n{N{k-l)/N{k))’ 


where eq{k) denotes the resulting error in the approximation of the control variable 
q and /c G {2, • ■ • ,8}. 
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DOES 

l|u - Uj7fc||L2(o) 

EOCn 

|7 - 75^JlL“(n) 

EOCy 

42 

0.0456202 

- 

0.3940558 

- 

146 

0.0259039 

-0.4542396 

0.1220998 

-0.9403796 

546 

0.0106388 

-0.6746618 

0.0356279 

-0.9338121 

2114 

0.0053128 

-0.5129453 

0.0104755 

-0.9042427 

8322 

0.0026798 

-0.4994327 

0.0030256 

-0.9063059 

33026 

0.0013372 

-0.5043272 

0.0008921 

-0.8860222 

131586 

0.0006675 

-0.5025385 

0.0002586 

-0.8957802 

525314 

0.0003340 

-0.5000704 

7.359881e-05 

-0.9077666 


Table 6.1 


Experimental order of convergence of scheme (l2.1QII - (l2.12l t when is used to approximate the 
solution to the optimization problem of Section p| with one observation point. The EOCq is in 
agreement with estimate (I4.13l t of Theorem \4-3\ the family {^k} quasi-uniform and, thus, h^^ ^ 
, which is what we observe. The EOCy reveals a quadratic order and illustrates that our 
error estimate II4.15II might be pessimistic. 


Table [6?n shows that, when approximating the optimal control variable, the EOCo 
is in agreement with the estimate (|4.13|1 . This illustrates the sharpness of the derived 
estimate up to a logarithmic term. We comment that, since the family is 

quasi-uniform, we then have that Consequently, (|4.13p reads as 

follows: 


l|u - U 5 .J|i 2 ( 0 ) < iV(/c) 2 |logiV(fc)|. (6.3) 

Table lOI also presents the EOCy obtained for the approximation of the optimal state 
variable y: k, N{k)~^. These results show that the derived error estimate (I4.15|) 

might be pessimistic. 

6.3. Optimization with point observations: four points. The objective 
of this numerical experiment is to test the performance of the fully discrete scheme 
(I2.10I1 - (I2.12I1 when more observation points are considered. 

Let us consider n = 2 and O = (0,1)^- The control bounds defining the set 
are given by a = —1.2 and b = —0.7. The state equation (14. 3p is replaced by (16.21) . 
which allows the incorporation of a forcing term f. We set 

Z = {(0.75,0.75), (0.75,0.25), (0.25,0.75), (0.25,0.25)}, 

with corresponding desired values 

7(0.75,0.75) = 1) 7(0.25,0.25) = 1 ) 7(0.75,0.25) = 7(0.25,0.75) = 0 - 5 - 

The exact optimal state variable is then given by 

y{xi,X2) = 2.75 — 2xi — 2x2 + 4x1X2 

and the exact optimal adjoint state is given by (EU). 

Table 16.21 shows that the EOCq is in agreement with estimate (14.1311 of Theo¬ 
rem Mi This illustrates the robustness of scheme (I2.10l) ~ (l2.12p when more observa¬ 
tions points are considered. 
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DOFs 

l|u - Uj7fc||L2(o) 

FOC 

|y - y^k\\L°°(n) 

FOC 

42 

0.0285416 

- 

0.0595256 

- 

146 

0.0285084 

-0.0009357 

0.0152388 

-1.0936039 

546 

0.0208153 

-0.2384441 

0.0039226 

-1.0288683 

2114 

0.0116163 

-0.4308717 

0.0010313 

-0.9868631 

8322 

0.0061821 

-0.4602926 

0.0002708 

-0.9758262 

33026 

0.0030792 

-0.5056447 

7.057710e-05 

-0.9755383 

131586 

0.0014908 

-0.5247299 

1.729492e-05 

-1.0173090 

525314 

0.0007618 

-0.4849766 

4.503108e-06 

-0.9720511 


Table 6.2 


Experimental order of convergence of scheme (l2.1QII - (l2.12l t when is used to approximate the 
solution of the problem of Section\^with four observation points. The EOCo is in agreement with 
estimate II4.13I I of Theorem \4.3\ the family {^k} quasi-uniform, so that h^^ ^ which 

is what we observe. The EOCy reveals a quadratic order and illustrates that our error estimate 
(I4.15l t might be pessimistic. 


6.4. Optimization with point observations: a three dimensional exam- 

pie. We set n = 3 and n = (0,1)^. We define a = —15 and b = —5. The optimal 
state is 


y(a;i,a;2,a:3) = ^^a;ia;2a;3(l - a;i)(l - X2)(l - xs), 
whereas the optimal adjoint state is defined by (EU. The set of observation points is 

Z = {(0.25,0.25,0.25), (0.25,0.25,0.75), (0.25,0.75,0.25), (0.25,0.75,0.75), 

(0.75,0.25,0.25), (0.75,0.25,0.75), (0.75,0.75,0.25), (0.75,0.75,0.25)} 

and we set yz = 1 for all z € Z. 


DOFs 

l|u - Uj7J|l2(q) 

FOCd 

1419 

0.0274726 

- 

3694 

0.0199406 

-0.3349167 

9976 

0.0120137 

-0.5100352 

27800 

0.0088690 

-0.2961201 

79645 

0.0067903 

-0.2537367 

234683 

0.0049961 

-0.2839348 

704774 

0.0036908 

-0.2753755 

2155291 

0.0026540 

-0.2950207 


Table 6.3 


Experimental order of convergence of scheme Il2.10ll - ll2.12ll when used to approximate the so¬ 
lution to the optimization problem of Section p| a three dimensional example. The EOCo is 
in agreement with estimate 114.1411 of Theorem \4.3\ the family is quasi-uniform and then 


Table 16.31 illustrates that EOCq is in agreement with estimate (14.141) of The¬ 
orem 231 due to the fact that the family {15^} is quasi-uniform and then Ri 

DOFs(fc)-i/3 

6.5. Optimization with singular sources. We now explore the performance 
of scheme (I2.10|) - (I2.12I) when is used to solve the optimization problem with singular 
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sources. We set n = 2 and ft = (0,1)^. We consider V = (0.5, 0.5) and the control 
bounds that define the set U 5 are a = 0.3 and b = 0.7. The desired state and the 
exact adjoint state correspond to 

p(a;i, a; 2 ) = —32a;ia:2(l — a;i)(l — X 2 ), Yd = — sin(27ra;) cos(27ra;) 

The exact optimal state is given by (EDI¬ 


DOFs 

l|u - Uj3?-J|L2(n) 

EOC 

30 

0.0940682 

- 

86 

0.0536485 

-0.5332256 

294 

0.0207101 

-0.7743303 

1094 

0.0068950 

-0.8369949 

4230 

0.0021408 

-0.8648701 

16646 

0.0006380 

-0.8836678 

66054 

0.0001850 

-0.8981934 

263174 

5.259841e-05 

-0.9098104 

1050630 

1.472536e-05 

-0.9196613 


TABLE 6.4 


Experimental order of convergence of scheme (l2.1Q|l - ()2.12ll when used to approximate the so¬ 
lution to the optimization problem with point sources of Section[^ The EOCo reveals a quadratic 
order and illustrates our error estimate (I5.2lt . 
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